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Abstract. 

We theoretically investigate the effect of intcrtube tunneling in topological 
superfluid phases of a quasi-one-dimensional Fermi gas with a Rashba-type spin- 
orbit interaction. It is shown that the effective Hamiltonian is analogous to that 
*^ , of a nanowire topological superconductor with multibands. Using a hidden mirror 

Ph ■ symmetry in the system, we introduce a new topological number that ensures the 

existence of non-Abelian Majorana zero modes even in the presence of intertubc 
tunneling. It is demonstrated from the full numerical calculation of self-consistent 
equations that some of Majorana modes survive against the intertube tunneling, when 
the number of one-dimensional tubes is odd in the y-direction. We also discuss a 
generalization of our consideration to nanowire topological superconductors. 

CO 
CN . 

O ■ 1- Introduction 

CO t 

Majorana fermions are real fermions which are equivalent to their own anti-particles. 
Since the pioneering works by Read and Green [1] and Kitaev [2], search for their 
elusive fermions opened an exciting new chapter in condensed matter physics. Original 
works in Refs. [1] and j2] predicted that the mysterious fermions exist as zero- 
energy quasiparticles bound at vortices and edges of a spinless p-wave superconductor. 
Subsequently, tremendous progress has succeeded in extending the platform for realizing 
Majorana fermions to some categories of so-called topological superconductors [U HJ |5l [6] . 
The remarkable consequence of the self-charge conjugate property of Majorana fermion 
is non-Abelian braiding statistics, where a pair of Majorana zero modes are created or 
annihilated by braiding their host vortices [7]. Hence, Majorana fermions possessing 
non-Abelian braiding statistics can provide a promising platform for fault-torelant 
topological quantum computation [8]. Moreover, it has recently been unveiled that zero- 
energy quasiparticles exhibit multifaceted properties, not only as a Majorana fermion 
but also as odd-frequency Cooper pair correlation [9j [10J [HJ [12} [13j [14] . 
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An ideal candidate for realizing non-Abelian Majorana zero modes was a chiral 
p-wave superconductor with half-quantum vortices, where the low-lying quasiparticles 
are effectively spinless [7J. A half-quantum fluxoid has been observed in mesoscopic 
annular rings of Sr 2 Ru0 4 [15], while the half-quantum vortices in superconductors 
are energetically unstable against integer vortices because of the absence of screening 
mechanism of spin current [16]. The another candidate of a chiral p-wave superfluid 
is the A-phase of superfluid 3 He confined to a restricted geometry with sub-micron 
thickness [17]. The thermodynamic stability of half-quantum vortices in such a 
superfluid is not trivial, because the Fermi liquid corrections which favor vortices with 
spin flow rather than mass flow are competitive to the strong coupling correction due 
to the spin fluctuation which stabilizes integer vortices without spin flow [IS] . 

On the contrary, it was demonstrated that a conventional s-wave superconductor 
can harbor non-Abelian Majorana zero modes [19j EUJ ETJ [22J E3]- The key finding 
lies in the two-dimensional Rashba-type spin-orbit interaction in background normal 
fermions, where the non-Abelian anyons are due to the phase twist of the spin- 
orbit interaction. The electron bands split by the spin-orbit interaction effectively 
convert the s-wave pairing to p ± ip pairing. A strong Zeeman field drives the 
quantum phase transition from a non-topological phase without non-Abelian anyons 
to a topological phase. This finding provides an another approach to the realization 
of topological quantum computation in condensed matters. Indeed, it has been 
proposed that topological superconductivity and Majorana fermions can be realized 
in a one- dimensional (ID) semiconducting wire proximity-coupled with an s-wave 
superconductor [2U I22J I2S1 |27J [281 I2H], where semiconductors, such as InSb, have 
very large g-factor and strong spin-orbit interaction. The signature of Majorana 
fermions has recently been observed through zero bias conductance peaks in a nanowire 
topological superconductor [301 Ell E2] and unconventional Josephson effect in hybrid 
superconductor-topological insulator devices [33] . 

Apart from superconducting materials, cold atoms with a p-wave Fershbach 
resonance [3^1 135] or with a synthetic gauge field [211 122] offer an alternative playground 
for Majorana fermions. A spin-orbit coupling with equal Rashba and Dresselhaus 
strengths can be synthetically induced by applying Raman lasers to atomic gases with 
hyperfine spin degrees of freedom [36], which has been implemented using 6 Li [37] 
and 40 K atoms [38] • In addition, schemes for creating Rashba and Dresselhaus spin- 
orbit coupling and three-dimensional analogue to Rashba spin-orbit coupling have 
theoretically been proposed [391 SD1 SI]- A ID geometry can be implemented by 
using a two-dimensional optical lattice, which was already utilized to search the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state in a spin-imbalanced Fermi gas [12]. Hence, 
cold atoms with Raman laser-induced spin-orbit coupling provide not only a promising 
platform for realizing Majorana fermions [211 EH1 HH 051 H6] but also an ideal system to 
study competition of various exotic superfluid phases including topological and FFLO 
phases. 

In this paper, we study topology and quasiparticle spectra of a quasi-lD Fermi gas 
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with a Rashba-type spin-orbit coupling. It has been demonstrated in Refs. [I3J 0U 0S] 
that a pure ID Fermi gas with spin-orbit coupling is accompanied by exactly zero-energy 
states bound at the end points of atomic clouds. However, since an actual experiment 
is performed on a bundle of weakly coupled tubes, intertube tunneling effects are not 
negligible. We here demonstrate that the intertube coupling plays an important role on 
determining the topological properties of Fermi gases and the effective Hamiltonian is 
analogous to nanowire topological superconductors with multibands [HI H3 SHI EH |5T] . 

This paper is organized as follows. We begin in Sec. 2 by introducing a tight- 
binding model for a bundle of ID Fermi gases with spin-orbit coupling. In Sec. 3, 
we clarify the topology of such a system, where a one-dimensional winding number is 
protected by a hidden mirror symmetry. It turns out that this system provides a cold 
atom analogue to nanowire topological superconductors with multibands. In Sec. 4, 
based on fully numerical calculations of self-consistent equations, we study the intertube 
tunneling effect on quasiparticle spectra. The final section is devoted to conclusions and 
discussions. Throughout this paper, we set h — k& — 1 and the repeated Greek indices 
imply the sum over x, y, z. 

2. Array of one-dimensional tubes 

We here consider a system composed of spin-orbit coupled two-component fermionic 
atoms with an s-wave attractive interaction, g. This is described by the Bogoliubov-de 
Gennes (BdG) Hamiltonian 



where denotes the Pauli matrices in spin space and the single-particle Hamiltonian 
density is defined as e(r) = — ;^V 2 — fJ> + V pot (r) — h^a^ with a confinement potential 
Vp Qt . The Zeeman field h = in Eq. (CQ) is naturally induced by implementing the 
spin-orbit coupling through two-photon Raman process [36] • In Eq. ([T|), S describes the 
spin-orbit coupling term, which is expressed in general as 



The quasiparticle states with the energy E n are obtained as the energy eigenstates 
of H(r), the so-called BdG equation, 1-L{r)ip n {r) = E n <p n (r), where <p = 
[u n ^, u n ^, i> n ,t, v n,i\ T denotes the wavefunction of quasiparticles. 

In order to isolate Majorana zero modes from the higher energy quasiparticle states, 
the fermionic atoms are confined by a two-dimensional optical lattice in the y-z plane 
in addition to a shallow harmonic potential along the x-direction, as shown in Fig. [TJ 
The system under this confinement potential is regarded as a two-dimensional array 
of N y x iV 2 -th ID tubes. The quasiparticle wavefunction is expanded in terms of the 
functions f ty {y) and / 4 (z) localized at (x,yt y ,z iz ) as (p(r) = Y,e.Vi{.x)fz y (y)fiA z ) 
with £ = (£ y ,£ z ). Employing the tight-banding approximation in the y-z plane, the 




(1) 



S(r) = iA^o-yd^. 



(2) 
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Figure 1. (a) Two-dimensional optical lattice potential and (b) schematic picture of 
the calculated system. 



Hamiltonian in Eq. ([T]) reduces to 



,A}(x) -§*_{x)-S*(x 



t£,£'{x) = ( -^-71 - (i + Vt(x) - h^a^ 8 tt t' 



•y 

where the e(x) is a 2N y N z x 2N y N z matrix given by 

2m dx 2 

— ty($£,£'+e y + ^£,£'-e v ) ~ tz($£,£'+e z + $£,£'-&*)■ (4) 

The effective potential Vi(x) is given as Vt(x) = \mu) 2 x 2 + V^. Here we have introduced 
£ = (^,4), e y = (1,0) and e z = (0, 1), where £ y = l,---,N y and £ Z = 1,---,N Z . The 
hopping energies between intertubes are denoted by t y and t z . The spin-orbit coupling 
term also reduces to 



St Ax) = ia 



+A Zfl (5e t £> + e z - A ztl 5e^-e z ) , (5) 

where A^ u describes an effective non-Abelian gauge field, A Vil = J f^(z)f^ +1 (y)A y ^fi y (y)fe z (z)dydz. 
Note that a quasi- ID Fermi gas has been reported in Ref. jl2] with a two-dimensional 
optical lattice potential. The number of tubes is typically about N y xN z ~ O(10 x 10), 
so the system should be treated finite system. 

Under the tight-binding approximation, the resulting BdG equation reduces to an 
effectively ID equation along the x-axis, 

KfA^tPnt'ix) = E n cp n/ (x), (6) 

which is numerically solved with the finite element method implemented with the 
discrete variable representation [5?]. The BdG equation (jHJ) is selfconsistently coupled 
with the gap equation for the pair potential, A^i/(x)=Ai(x)5^i/, 

A A%) =dYl [ U n,tA X ) V n,lA X )f( E ") + U nAA X ) V hA X )f(~ E n)] » ( 7 ) 

where f(E) = l/(e E / T + 1) is the Fermi distribution function at a temperature T. In 
addition, the chemical potential fi is determined so as to preserve the total particle 
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number 

N=J2 fdr [\u n , a , £ (x)f f(E n ) + \v n , a>t (x)\ 2 f(-E n )] . (8) 

The sum in Eqs. ([7]) and ([8]) is taken over E n G [0,^], where E c denotes the energy 
cutoff. The effective ID coupling constant, g, in Eq. ((Tj) is expressed in terms of an 
effective ID scattering length aio as [52] 

2 2 a 3D , x 

9 = = -j t 9 

maiD m Aa 3 D — d± 

where a 3 D denotes a three-dimensional scattering length and d± is the characteristic 
harmonic oscillator length in y- and z-axes. The ID coupling constant g is always 
negative for a negative a 3 D, such as 6 Li atoms. The constant A is given as A ^1.0326. 
It is seen from Eq. Q that the pairing interaction, g, can be controlled by changing the 
characteristic length scale of the confinement potential in the y-z plane, d±. 

3. Topology of the effective Hamiltonian 

We here consider a two-dimensional Rashba-type spin-orbit interaction, 

A^Oyd^ = K x a y d x - K y o x d y . (10) 

To clarify the topological property of the effective Hamiltonian, we here ignore the 
shallow trap potential along x-direction, i.e., u x — > 0. Then, the effective Hamiltonian 
in Eq. ([3]) is rewritten with —id x — >• k as 



^Tt\k) = ef^,{k) - {h x a x + h z a z + n x ka y } 8 tjl 

- Ai^ayTy + \hy<T y 8ifi + ik y (T x [d ^(j +&y - S^^}] T Q , (11) 

where ef\, describes the single-particle Hamiltonian density without the Zeeman term 
and r M denotes the Pauli matrices in particle-hole space. Here, without the loss of 
generality, is assumed to be real. We have also assumed h y = 0. It is seen 
from Eq. (TTTT) that the array of ID tubes with a spin-orbit interaction is analogous 
to a semiconductor-superconductor nanowire with iV-th electron bands [51], where 
N = N y xN z . 

We find that our system supports two different topological numbers. The first one 
is the ID Z 2 topological number for class D topological phases: Because of superfluidity, 
the BdG Hamiltonian (fTT|) has the particle-hole symmetry, which allows us to define the 
ID 7*2 topological number. The ID Z 2 number is defined as 

i r 

- / dkA{k) +mod. 2, (12) 
where A(k) is the geometrical phase, 

E n {k)<0 t 



V 
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with \(p n ^(k)) the Bloch wave function of an negative energy state of "H eff (/c). When v 
is odd (even), the system is topologically non-trivial (trivial). 

The second topological number comes from a remnant of a mirror reflection 
symmetry of the system. If one temporary neglects the Zeeman fields h, our system is 
invariant under the mirror reflection to the zx-plane, as well as the time-reversal. Once 
the Zeeman fields are applied, the mirror symmetry is lost, but a combination of the 
mirror reflection and the time- reversal is still preserved if h y — 0. Consequently, the 
Hamiltonian H e s(k) with h y = holds the following Z<i symmetry, 

rM zx n cS (k)Mlr^ = n* eS (-k), (u) 

where T=io y K is the time-reversal operator with the complex conjugate operator K, 
and M. zx = lOyU the mirror reflection operator. Here U is the operator flipping the 
y-component of I = (£ y ,£ z ), 

Uz,i> = St v ,N y +i-z' y h z ,t' z - (15) 

Combining the Z 2 symmetry with the particle-hole symmetry, CH c fi(k)C~ l = 
—%* eS {—k), we define the chiral symmetry operator, Y = TM. zx C = t x U , and it turns out 
that T is anti-commutable with the effective Hamiltonian 

{T,n cS (k)} = 0. (16) 

This implies that the BdG Hamiltonian 1-L c s{k) holds the chiral symmetry. Then, the 
ID winding number is defined as [531 EH EHl [56] 



w 



1 f°° 

— J dktr [rn;i(k)d k n c{i (k)] , (17) 



47ri 

which takes an integer. 

As a consequence of the bulk-edge correspondence, these two ID topological 
numbers ensure the existence of zero energy states appearing in the end points of ID 
tubes. Here we note that the parities of these two topological numbers coincide with 
each other, 

(-ir=(-ir, as) 

which implies that w can be nonzero even when v is trivial, but the opposite is not true. 
Therefore, the actual number of the zero energy states is determined by w unless the 
Z2 symmetry (1141) is broken macroscopically. In addition, the particle-hole symmetry 
of superfluidity results in the Majorana property of the zero energy state, where the 
creation operator 7^ is equivalent to its own annihilation, 7^=0 = 7^ [511 EH]. In 
summary, the winding number w ensures the existence of Majorana zero modes, whereas 
once the Z2 symmetry is broken by turning on h y for example, w becomes ill-defined 
and the ID Z 2 number v in Eq. (1121) determines the topological stability of the Majorana 
zero modes. 

In the case of a two-dimensional disk geometry, where the two-dimensional optical 
lattice potential is absent, the topological property is characterized by the first Chern 
number [22] - The nontrivial value of the Chern number ensures the existence of the 
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w 



gapless state localized at the circumference of the atomic cloud. In addition, the edge 
states carry the macroscopic mass current. 

4. Intertube tunneling effect 

Let us start with a pure ID system with t y = t z = k y = 0. Here, the BdG equation 
OH]) coupled with the gap equation (JTJ) is numerically solved with the set of parameters: 
T = 0, E c = 4Ep, and k—1. The Zeeman field is applied along z-axis, which does not 
break the Z 2 symmetry in Eq.(IT4l): h = (0,0, h). In realistic situations [36j HO], the 
strength of the spin-orbit coupling, k x , depends on the wave length of applied lasers. 
Throughout this work, we fix the pairing interaction, 7 = 1.4, where 7 = ^^(^~) is the 
dimensionless coupling constant [52] , where d = ^Jl/{mu x ) is the harmonic oscillator 
length. The total particle number is fixed to be iV = 200 in each tube, where the Fermi 
energy per one tube is E-p = 100oj x . 

For the pure ID system, topology of the BdG Hamiltonian for each tube is 
characterized by the winding number w in Eq. ( ITT)) with N y = N z = 1. With spatially 
constant A and fi, w is given by 

f dkd k \n[detQ(k)}, (19) 

where Q(k) = e^°\k) — (h x a x + h z a z + nka y ) + iAa y . Then, it is found that w — 1 when 
\h\ > h c = a//i 2 + A 2 . Therefore the lower magnetic field regime is the topologically 
trivial phase and the critical field at h c involves the topological phase transition. For 
our system, however, a more careful consideration is needed. Since fermionic atoms 
are confined by a trap potential, a spatially inhomogeneous superfluid is realized 
naturally. In addition, the inhomogeneous pair potential A(x), which is self-consistently 
determined by the gap equation and the BdG equation, depends on the Zeeman fields 
significantly. In contrast to semiconductor-superconductor junction systems, these 
two characteristics cannot be neglected. This means that within the local density 
approximation, the critical field h c and ID winding number w should be replaced by 
h c {x) = a/ j2 2 (x) + A 2 (x) and w(x), where is the local chemical potential including 
the confinement potential along the x-axis. The inhomogeneity and self-consistency of 
A (a;) and /i(x) play a critical role on the topological property. 

Figure [2(a) shows the spatial profile of A(x) at h = 0.20Ep and 0.36-Ef, where the 
harmonic trap potential along the x-axis, ^mu^x 2 , is taken into account. In the case of 
h = 0.36E F , the intermediate region between 8.8cZ< \x\ <14.6d, where the pair potential 
has a dip, becomes topologically non-trivial, i. e. it satisfies h > h c (x) and w(x) 7^ 0, while 
the inner region within |x| <8.8c? is not. It is also found that the outermost region is not 
topological again, because the local chemical potential, /i(x)=fj, — jmu x x 2 , changes its 
sign at the Thomas-Fermi edge x = xtf = dy / 2fi/uj x ~ 14c? and its magnitude becomes 
large so as h < h c (x). 

In the absence of the hopping between tubes, four zero modes appear in each 
tube. At h — 0.3QEp, \w\ = 1 can be realized in the region within 8.8<i < \x\ < 14.6d, 
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Figure 2. (a) Spatial profiles of the pair potential A(x) in a pure ID system 
with t y = t z = and h y = 0. (b) Wavefunctions of the lowest (upper panel) and 
second lowest (lower panel) energy eigenstates at h = 0.36E-p, where the energies are 
£ = 3.77xlCT 10 i; F and 7.28 x 1Q~ 9 E F . 




0.4 0.45 0.5 



! ! I I I 



2 3 4 5 

N„ 



Figure 3. Field-dependence of quasiparticle energy spectra in a pure ID case (a) 
and quasi-lD cases (b). (c) Quasiparticle spectra as a function of N y at h — 0.36-Ep 
which corresponds to the topological phase in the case of N y = 1 . In all the data, the 
hopping t y and the strength of the spin-orbit interaction h y is set to be t y = 0.01-Ef 
and k y j ' k x = 0.5. 



otherwise w = 0. As shown in Fig. Mb), the lowest and second lowest energy states 
are tightly bound at the phase boundaries at x ~ ±8d and ±15<i. We summarize 
the field-dependence of the energy spectrum E n in Fig. [3]^a). It is seen that the low- 
lying eigenenergies go to zero and yield the oscillation as h increases, because of the 
interference between two zero modes localized at x/c?~±8 and ±15 [5T| [58]. 

Now, let us clarify how the intertube tunneling affect the low-lying quasiparticle 
spectra. The bundle of one-dimensional tubes are coupled with each other through the 
hoppings t y and t z and the spin-orbit interactions k y . In this situation, the topological 
winding number w is not defined for each tube, but is defined only for a whole system 
of tubes. As a result, some of Majorana zero modes become non zero modes as is 
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Figure 4. Local density of states J\f a ^(x,E) denned in Eq. (|20p in the case of N y = 3 
and h = 0.36-Ef- The LDOS at £ y = 1 for a =t an( i i is displayed in (a) and (b), 
respectively, (c) and (d) shows the LDOS for a =j" and \. at t y = 2. The other 
parameters are same as those in Fig. [3] 

shown below. Note that if the hopping t z is small enough, it does not change the 
winding number and the topological property because of the two- dimensionality of the 
Rashba-type spin-orbit coupling. This implies that zero energy states are dispersionless 
and insensitive to N z . Figures EJ^b) and (c) show the field-dependence of low-lying 
quasiparticle energies for N y = 2 and 3, where we fix the parameters, t y = O.OIE'f and 
k v /k x = 0.5. In the case of N y — 2, the low-lying energies are sensitive to the intertube 
tunneling effect, t y and k y , where the energies of the quasiparticles bound at the phase 
boundaries increase as t y and k y increase. In contrast, as shown in Fig. [3^b), the low- 
lying quasiparticle states in the case of N y = 3 is unchanged by the intertube tunneling 
effect and become the zero energy as h increases. Figure [3]^c) summarizes the low-lying 
spectra as a function of N y , where h is fixed to be hj Ep = 0.36 which corresponds to the 
topological phase in the case of N y — 1. It is clear that the zero energy states exist when 
N y is odd, while in the case of even A^'s, the intertube tunneling lifts the low-lying 
states from the zero energy. 

To understand the low-lying energy states in the case of odd N y s, which remain 
zero modes, we show in Fig. H] the local density of states (LDOS) at the site £ defined 
as 

AT^(x,E)= [\u n ^{x)\ 2 5{E-E n ) + \v n ^{x)\ 2 5{E + E n )]. (20) 

E n >0 
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Figure 5. Amplitude of wave functions \<PE=o,e(x)\ for the lowest energy state in the 
case of (a)N y = 3, (b)N y = 5, and (c)N y = 7, respectively. The parameters are same 
as those in Fig. [3] 

The LDOS for up-spins (down-spins) in the case of N y = 3 is displayed in Figs. Ufa) and 
(c) ((b) and (d)). It is clearly seen from Figs. IU(a) and (b) that the LDOS for up-spins 
at i y — 1 and 3 is accompanied by the zero energy states which are bound at the end 
points of the tube [xjd ~ ±15), while A/4, has sharp peaks at x/gZ~±8 corresponding to 
the inner phase boundaries between w(x) —0 (non-topological region) and 1 (topological 
region) in the case of N y — 1. In contrast, the intertube tunneling split the zero energy 
states to the positive and negative energy states in the LDOS at £ y = 2, which have the 
mini-gap with ±0.02Ep. In Figj5], we also show the amplitude of the zero modes for 
N y — 3, Ny — 5 and N y — 5, respectively. In all cases, the wave functions have large 
amplitudes at the tubes located at odd £ y s. At even £ y s, their amplitudes are almost 
negligible. 

These intertube tunneling effects can be understood as follows. As was shown 
above, when one neglects the intertube couplings t y , t z and k y , each tube supports four 
Majorana zero modes localized at x ~ ±8d and x ~ ±15gL Now let us denote one of 
them (say, the zero mode localized at x ~ ±15g/) as (£ y — 1, ■ • • , N y ), and consider 
how the intertube couplings affect on them. When t y and k y are turned on, the zero 
modes on neighboring tubes are coupled by the intertube tunneling, 

% = ^7 l72 + zt7 2 73 H h itjNy-ilNy, (21) 

where t denotes the induced tunneling coupling. Note that t is real since 7^ is a 
Majorana zero mode satisfying = 7] . Equation ( 12TT) is rewritten as H = 
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with 



/ 



H 



\ 



-it 
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-it 
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it 




(22) 



/ 



and T = (71, 72, • • • , jN y -i, lN y Y- Diagonalizing the N y x N y matrix H, one can examine 
the effects of the intertube tunneling. 

It can be easily shown that H has a single zero eigenvalue for odd iV y 's, while it 
does not have for even N y 's. This result naturally explain why Majorana zero modes 
survive only for odd N y 's. One also finds that the zero eigenstate of H has the following 
form 



[1,0 

;i,o 

'1,0 



for N„ 



1,0, l) 4 , 
1,0,1,0, 



for N„ 



)*, for N y = 7, (23) 

which explains qualitatively why the remaining Majorana zero modes illustrated in Figj5] 
have large amplitudes on tubes at odd £ y s. 

The robustness of the zero-energy states against the intertube tunneling is also 
understood by the topological number w. As we mentioned above, even in the presence 
of intertube tunneling, the winding number w is well-defined for a whole system of tubes. 
Since one can turn off the intertube tunneling without the bulk gap closing, the value 
of w can be evaluated by setting t y = t z = k y = in Eq. ffTTl) . Then one obtains 



w 



1 f°° 

trU / alkd k In [det Q(k)] . 

2m 



Noting 



tvU = 
one can valuate w as 



\w\ 



for even N y 's 
N z for odd Ays ' 



for odd ^'s 
0, for even iV^'s 



(24) 



(25) 



(26) 



when \/hl + h 2 z > h c . This implies that Majorana zero modes survive for odd A^'s. 



5. Conclusions and Discussions 



In this article, we have investigated the effect of intertube tunneling in a quasi-one- 
dimensional Fermi gas with a Rashba-type spin-orbit coupling. From the argument 
based on the symmetry of the effective Hamiltonian, the topological property has been 
studied. In the absence of the intertube tunneling, the one- dimensional winding number 
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in Eq. (j!7p ensures the existence of Majorana zero modes bound at the end points of 
each tube. It also ensures the existence of the Majorana zero modes in the presence 
of the intertube tunneling if the number of tubes in the ^-direction is odd. Using 
full numerical calculations of self-consistent equations, we have confirmed that this 
topological property is clearly reflected in low-lying quasiparticle states. These behaviors 
of low-lying quasiparticles associated by the Z 2 -symmetry protected topology might be 
detectable through the momentum- resolved radio-frequency spectroscopy [59], [60] . 

Finally, we would like to mention about a generalization of the present consideration 
to semiconductor-superconductor nanowire systems. The one-dimensional winding 
number ( fTTl) introduced in this paper is also applicable to semiconductor-superconductor 
nanowire with multichannels. If we consider the nanowire extending in the ^-direction on 
the top of an s-wave superconductor in the x?/-plane, the system is naturally supposed 
to be invariant under the mirror reflection, y — > —y, to the xz-plane. This mirror 
symmetry could be broken under Zeeman fields, but the Z 2 symmetry (fh4"j) remains if 
the Zeeman fields are applied in the x- or z-direction. Then, the topological number 
w in Eq. (fTTl) is defined in the same manner. From arguments similar to the above, 
one finds that w is nonzero if the Zeeman filed h satisfies \h\ > h c and the number of 
the channels in the y-direction of the nanowire is odd. Indeed, under this condition, 
\w\ is equal to the number of channels in the ^-direction of the nanowire. Note that, in 
contrast to the ID Z 2 number in Eq. (fl~2l . w can be nonzero even when the total number 
of channels in the nanowire is even, since it is given by sum of the channels in the z 
and ^/-directions. As well as the fermionic gas system studied in this paper, the local 
density operator of the Majorana zero modes vanishes [55], and thus the Majorana zero 
are stable against weak non-magnetic disorders. 

Acknowledgments 

The authors are grateful to Takuto Kawakami for fruitful discussions and comments. 
This work was supported by a Grant-in-Aid for Scientific Research from MEXT of 
Japan, "Topological Quantum Phenomena" No. 22103005 and JPSJ No. 22540383. 

References 

[1] N. Read and D. Green, Phys. Rev. B 61 (2000) 10267. 
[2] A. Y. Kitaev, Phys. Usp. 44 (2001) 131. 

[3] Y. Tanaka, M. Sato, and N. Nagaosa, J. Phys. Soc. Jpn. 81 (2012) 011013. 
[4] X.L. Qi and S.C. Zhang, Rev. Mod. Phys. 83 (2011) 1057. 

[5] A.P. Schnyder, S. Ryu, A. Furusaki, and A.W.W. Ludwig, Phys. Rev. B 78 (2008) 195125. 
[6] A. Kitaev, AIP Conf. Proc. 1134 (2009) 22. 
[7] D. A. Ivanov, Phys. Rev. Lett. 86 (2001) 268. 

[8] C. Nayak, S.H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Rev. Mod. Phys. 80 (2008) 
1083. 

[9] T. Yokoyama, Y. Tanaka, and A. A. Golubov, Phys. Rev. B 78 (2008) 012508. 
[10] Y. Tanuma, N. Hayashi, Y. Tanaka, and A.A. Golubov, Phys. Rev. Lett. 102 (2009) 117003. 



Topological phases of quasi- one- dimensional fermionic atoms 



13 



[11] S. Higashitani, S. Matsuo, Y. Nagato, K. Nagai, S. Murakawa, R. Nomura, and Y. Okuda, Phys. 
Rev. B 85 (2012) 024524. 
T. Daino, M. Ichioka, T. Mizushima, and Y. Tanaka, Phys. Rev. B 86 (2012) 064512. 
Y. Tsutsumi and K. Machida, J. Phys. Soc. Jpn. 81 (2012) 074607. 
Y. Asano and Y. Tanaka, larXiv: 1204.42261 

J. Jang, D.G. Ferguson, V. Vakaryuk, R. Budakian, S.B. Chung, P.M. Goldbart,and Y. Maeno, 

Science 331 (2011) 186. 
S.-K. Chung, H. Bluhm, and E.-A. Kim, Phys. Rev. Lett. 99 (2007) 197002. 
Y. Tsutsumi, T. Mizushima, M. Ichioka, and K. Machida, J. Phys. Soc. Jpn. 79, (2010) 113601. 
T. Kawakami, T. Mizushima, and K. Machida, J. Phys. Soc. Jpn. 80 (2011) 044603. 
M. Sato, Phys. Lett. B 575 (2003) 126. 
L. Fu and C.L. Kane, Phys. Rev. Lett. 100 (2008) 096407. 
M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. Lett. 103 (2009) 020401. 
M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. B 82 (2010) 134521. 

C. Zhang, S. Tewari, R.M. Lutchyn, and S. Das Sarma, Phys. Rev. Lett. 101 (2008) 160401. 
J.D. Sau, R.M. Lutchyn, S. Tewari, and S. Das Sarma, Phys. Rev. Lett. 104 (2010) 040502. 
R.M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. Lett. 105 (2010) 077001. 
J. Alicea, Phys. Rev. B 81 (2010) 125318. 

Y. Orcg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 105 (2010) 177002. 
P.W. Brouwer, M. Duckheim, A. Romito, and F. von Oppen, Phys. Rev. B 84 (2011) 144526. 
J. Alicea, Y. Orcg, G. Refael, F. von Oppen, and M. P. A. Fisher, Nat. Phys. 7 (2011) 412. 
V. Mourik, K. Zuo, S.M. Frolov, S.R. Plissard, E.P.A.M. Bakkers, and L.P. Kouwenhoven, Science 
336 (2012) 1003. 

M.T. Deng, C.L. Yu, G.Y. Huang, M. Larsson, P. Caro, and H.Q. Xu, Nano Lett. 12 (2012) 6414. 

A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum, and H. Shtrikman, Nat. Phys. 8 (2012) 887. 
J.R. Williams, A.J. Bestwick, P. Gallagher, S.S. Hong, Y. Cui, A.S. Blcich, J.G. Analytis, I.R. 

Fisher, and D. Goldhaber-Gordon, Phys. Rev. Lett. 109 (2012) 056803. 
T. Mizushima, M. Ichioka, and K. Machida, Phys. Rev. Lett. 101 (2008) 150409. 
V. Guraric and L. Radzihovsky, Ann. Phys. (NY) 322 (2007) 2. 
Y.-J. Lin, K. Jimcnez-Garci'a, and LB. Spiclman, Nature (London) 471 (2011) 83. 
L.W. Chcuk, A.T. Sommcr, Z. Hadzibabic, T. Ycfsah, W.S. Bakr, and M.W. Zwierlcin, Phys. Rev. 

Lett. 109 (2012) 095302. 
P. Wang, Z.-Q. Yu, Z. Fu, J. Miao, L. Huang, S. Chai, H. Zhai, and J. Zhang, Phys. Rev. Lett. 
109 (2012) 095301. 

G. Juzeliunas, J. Ruscckas, and J. Dalibard, Phys. Rev. A 81 (2010) 053403. 

D. L. Campbell, G. Juzeliunas, and I. B. Spiclman, Phys. Rev. A 84 (2011) 025602. 

B. M. Anderson, G. Juzeliunas, V.M. Galitski, and I. B. Spielman, Phys. Rev. Lett. 108 (2012) 
235301. 

Y.A. Liao, A.S.C. Rittncr, T. Paprotta, W. Li, G.B. Partridge, R.G. Hulet, S.K. Baur, and E.J. 

Mueller, Nature (London) 467 (2010) 567. 
X.-J. Liu and H. Hu, Phys. Rev. A 85 (2012) 033622. 
X.-J. Liu and P.D. Drummond, Phys. Rev. A 86 (2012) 035602. 
R. Wei and E.J. Mueller, Phys. Rev. A 86 (2012) 063604. 

H. Zhai, Int. J. Mod. Phys. B 26 (2012) 1230001. 
A.C. Potter and P.A. Lee, Phys. Rev. Lett. 105 (2010) 227003. 

R. M. Lutchyn, T.D. Stanescu, and S. Das Sarma, Phys. Rev. Lett. 106 (2011) 127001. 
T.D. Stanescu, R.M. Lutchyn, and S. Das Sarma, Phys. Rev. B 84 (2011) 144522. 
G. Kclls, D. Meidan, and P. W. Brouwer, Phys. Rev. B 85 (2012) 060507(R). 
S. Tewari, T.D. Stanescu, J.D. Sau, and S. Das Sarma, Phys. Rev. B 86 (2012) 024504. 
X.-J. Liu, H. Hu, and P.D. Drummond, Phys. Rev. A 76 (2007) 043605. 
M. Sato and S. Fujimoto, Phys. Rev. B 79 (2009) 094504. 



Topological phases of quasi- one- dimensional fermionic atoms 



11 



[54] M. Sato, Y. Tanaka, K. Yada, and T. Yokoyama, Phys. Rev. B 83 (2011) 224511. 

[55] T. Mizushima, M. Sato, and K. Machida, Phys. Rev. Lett. 109 (2012) 165301. 

[56] X.-G. Wen and A. Zee, Nucl. Phys. B316 (1989) 641. 

[57] T. Mizushima and K. Machida, Phys. Rev. A 82 (2010) 023624. 

[58] M. Cheng, R.M. Lutchyn, V. Galitski, and S. Das Sarma, Phys. Rev. Lett. 103 (2009) 107001. 

[59] J.T. Stewart, J.P. Gaebler, and D.S. Jin, Nature (London) 454 (2008) 744. 

[60] S.-G. Peng, X.-J. Liu, H. Hu, and K. Jiang. larXiv:1210.2160l 



